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ABSTRACT 



Over the last decade. X-ray observations unveiled the existence of several classes of isolated neutron stars (INSs) which are radio- 
quiet or exhibit radio emission with properties much at variance with those of ordinary radio pulsars. The identification of new 
sources is crucial in order to understand the relations among the different classes and to compare observational constraints with 
theoretical expectations. A recent analysis of the 2XMMp catalogue provided less than 30 new thermally emitting INS candidates. 
Among these, the source 2XMM J104608.7-594306 appears particularly interesting because of the softness of its X-ray spectrum, 
kT = 117 ± 14eV and A'h = (3.5 ± 1.1) x 10^' cm"^ (3(t), and of the present upper limits in the optical, iub i 26, my ~ 25.5 and 
Tor ^ 25 (98.76% confidence level), which imply a logarithmic X-ray-to-optical flux ratio log(Fx/Fv) S 3. 1 , corrected for absorption. 
We present the X-ray and optical properties of 2XMM J104608. 7-594306 and discuss its nature in the light of two possible scenarios 
invoked to explain the X-ray thermal emission from INSs: the release of residual heat in a cooling neutron star, as in the seven radio- 
quiet i?05Ar-discovered INSs, and accretion from the interstellar medium. We find that the present observational picture of 2XMM 
J104608. 7-594306 is consistent with a distant cooling INS with properties in agreement with the most up-to-date expectations of 
population synthesis models: it is fainter, hotter and more absorbed than the seven ROSAT sources and possibly located in the Carina 
Nebula, a region likely to harbour unidentified cooling neutron stars. The accretion scenario, although not entirely ruled out by 
observations, would require a very slow (~ 10km s"') INS accreting at the Bondi-Hoyle rate. 

Key words, stars: neutron ~ X-rays: individuals: 2XMM 1104608.7-594306 - Catalogs 



1. Introduction 

A major outcome of the ROSAT mission has been the discov- 
ery of a group of radio-quiet isolated neutron stars (INSs) with 
properties clearly at variance from the bulk of the population 
of rotation-powered pulsars. At present, the group comprises 
seven objects (whence the nickname "Magnificent Seven", or 
M7) sharing similar properties, which include purely thermal X- 
ray spectra with very soft temperatures (kT ~ 40-lOOeV) and 
low absorption column densities (A^h ~ few 10^" cm"^), long 
spin periods (P ~ 3-10 s), intense magnetic fields (B ~ 10'^- 
10^"* G), very faint optical counterparts (wb ^ 25 when detected, 
implying X-ray-to-optical flux ratios Fx/Fopi ^ 10"*) and no as- 
soc iations with supernova remnants or stellar companions (see 
e.g. lHaberill2007l for a recent review). 

For some tim e, f ollowing the pio neering work by 
lOstriker et all ( Il970h and lShvartsmanI (1197 Ih . it was thought that 
these sources might belong to the long sought population of old 
INSs powered by accretion of the interstellar medium (ISM), 



Send offprint requests to: A. M. Pires 

* Based on the public data archive of XMM-Newton an ESA Science 
Mission with instruments and contributions directly funded by the 
ESA Member States and the USA (NASA) and a Chandra Legacy 
programme. Optical observations were performed at the European 
Southern Observatory, Paranal, Chile, under programme ID 079. D- 
0633(A), and at the Southern Astrophysical Research Telescope, Cerro 
Pachon, Chile. 



which was expected to be detected by ROSAT (see iTreves et al.l 
2000, for a re view). However, the me asurement of large proper 
motions (e.g. Mot ch et al.ll2007l '2008, and references therein) 
strongly supports t he idea that the M7 represent a local popu- 
lation (d < 500 pc. iPosselt et al.ll2007l Ivan Kerkwiik & KaplanI 
'2007') of middle-aged neutron stars (~ 10^-10* yr), which give 
off thermal radiation as they cool down. The similar spin periods 
and their position on the P-P diagram raised the possibility that 
they could be related to other populations of INSs - like magne- 
tars, high magnetic field radio pulsars (HBPSRs ) and rotating 
radio transients (RRATs; iMcLaughlin et al.ll2006h - which is an 
intrigui ng and still debated is sue (e.g. lPopov et a'Ll l2006; Kapla^ 
120081: iKeane & Krameij|2008l) . 

Considering that within 1 kpc the M7 appear in compara- 
ble numbers as y oung (^ few Myr) radio and y-ray pulsars 
jPopov et alj2003b . they may represent the only identified mem- 
bers of a large, yet undetected, elusive population of radio-quiet 
and thermally emitting INSs. The discovery of new candidates is 
then mandatory in order to make any progress towards an under- 
standing of their properties as a population and of their relations 
with other classes of Galactic INSs. 

Despite the relatively small field of view and low sky cov- 
erage of the XMM-Newton ObservatorjO], its large effective area 
and good positional accuracy at soft X-ray energies make it ideal 



' In a near future, eROSITA will become an excellent tool to search 
for INSs (http : //www . mpe . mpg . de/pro j ects . htinl#erosita). 
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OBSID 


Observation Date 


Detector 


Filter 


texp 


B 


teff 


Counts 


Rate 


NearG 










(ks) 


(arcmin) 


(ks) 




(10"- s"') 




112580601 


2000-07-26 


pn 


thick 


27.7 


9.7 


15.0 


558 


3.73 ±0.18 


yes 






Ml 


thick 


33.2 


8.7 


17.8 


244 


1.37 ±0.10 


no 






M2 


thick 


30.1 


9.1 


15.7 


198 


1.26 ±0.10 


yes 


112580701 


2000-07-27 


pn 


thick 


8.0 


9.7 


4.3 


161 


3.7 ±0.3 


yes 






Ml 


thick 


10.9 


8.7 


5.9 


88 


1.50 ±0.19 


no 






M2 


thick 


7.9 


9.1 


4.1 


41 


1.00 ±0.18 


yes 


112560101 


2001-06-25 


pn 


thick 


21.5 


15.8 


6.7 


291 


4.4 ±0.3 


yes 


112560201 


2001-06-28 


pn 


thick 


20.5 


15.9 


6.3 


200 


3.2 ±0.3 


yes 


145740101 


2003-01-25 


Ml 


thick 


6.9 


8.3 


4.6 


78 


1.68 ±0.22 


yes 


145740201 


2003-01-27 


Ml 


thick 


6.9 


8.3 


4.6 


73 


1.59 ±0.21 


yes 






M2 


thick 


6.9 


8.0 


4.8 


69 


1.44 ±0.20 


yes 


145740301 


2003-01-27 


Ml 


thick 


6.8 


8.3 


4.5 


59 


1.30 ±0.20 


yes 






M2 


thick 


6.8 


8.0 


4.7 


68 


1.44 ±0.19 


yes 


145740501 


2003-01-29 


M2 


thick 


6.9 


7.9 


4.8 


67 


1.41 ±0.19 


yes 


160160101 


2003-06-08 


Ml 


thick 


16.5 


8.5 


8.8 


127 


1.44 ±0.15 


yes 






M2 


thick 


15.5 


8.2 


8.6 


114 


1.33 ±0.15 


no 


160160901 


2003-06-13 


Ml 


thick 


31.1 


8.5 


16.6 


207 


1.25 ±0.10 


no 






M2 


thick 


31.1 


8.3 


17.1 


217 


1.27 ±0.10 


no 


145780101 


2003-07-22 


Ml 


thick 


8.4 


9.4 


4.6 


76 


1.66 ±0.21 


no 


160560101 


2003-08-02 


M2 


medium 


11.8 


9.1 


6.3 


121 


1.93 ±0.19 


no 


160560201 


2003-08-09 


Ml 


thick 


12.2 


8.6 


6.8 


105 


1.54 ±0.17 


no 






M2 


medium 


12.1 


9.1 


6.4 


96 


1.50 ±0.18 


no 


160560301 


2003-08-18 


Ml 


thick 


18.5 


8.6 


10.8 


155 


1.42 ±0.13 


no 






M2 


medium 


18.5 


9.4 


10.0 


160 


1.60 ±0.14 


no 


206010101 


2004-12-07 


pn 


medium 


19.3 


17.3 


5.8 


366 


6.3 ±0.4 


no 


311990101 


2006-01-31 


pn 


thick 


24.3 


7.7 


15.9 


837 


5.27 ±0.26 


no 






M2 


thick 


65.3 


8.2 


44.2 


731 


1.65 ±0.15 


no 


9488 


2008-09-05 


ACIS-I 


OBF 


60.0 


6.7 


57.0 


567 


0.95 ± 0.04 


no 



Table 1. Description of the set of XMM-Newton and Chandra observations that serendipitously detected source XMM J1046. The 
exposure times (texp), filtered for background flares, and the effective exposures (teff) on XMM J 1046, accounting for vignetting, 
are reported; is the source off-axis angle. Net source photons and the on-axis count rates are in the 0.15-3keV (XMM-Newton) 
and 0.5-3 keV (Chandra) energy bands. 



to look for faint INS candidates. Recentlv. lFires & MotchI (|2008|) 
reported on the preliminary results of a programme aimed at 
identifying new thermally emitting INSs in the 2XMMp cata- 
loguq3 This version of the catalogue, released on July 2006, 
contains more than 120 thousand sources, covering a total sky 
area of ~ 285 deg^. Out of the considered ~ 7.2x10"^ EPIC pn 
sources with count rates above 0.01 s"', less than 30 good can- 
didates met all the selection criteria - i.e. well detected point- 
like sources showing soft thermal spectrum (kT < 200 eV and 
A^H - 10'^ - 10^^ cm"^), with no optical candidates and no iden- 
tifications in over 170 astronomical catalogues. The brightest 
and most promising INS candidate is source 2XMM J104608.7- 
594306 (hereafter XMM J1046). In this work we discuss the 
properties of XMM J 1046 in connection with the two main op- 
tions concerning its nature, a cooling or an accreting INS. The 
description of the selection procedure and the results of follow- 
up optical observations on a sample of the X-ray brightest INS 
candidates will be presented elsewhere (Pires et al., in prepara- 
tion). 

2. XMM J1046 in X-rays and the optical 

2.1. Observations and data reduction 

Thanks to the fact that source XMM J 1046 is at an angular dis- 
tance of only ~ 8.5 arcmin from the well studied binary system 
Eta Carinae, it was serendipitously observed by the EPIC pn 



and MOS detectors on board XMM-Newton on many dififerent 
occasions. As part of an observational campaign to study the 
star-forming region of the Carina Nebula, XMM J 1046 was also 
serendipitously observed once in a recent Chandra observation 
(September 2008). Although this X-ray mission has observed 
Eta Carinae many times, XMM J 1046 was not detected before 
given the smaller field of view of the ACIS instruments relative 
to the EPIC cameras. 

We analysed 16 XMM-Newton archival observations, span- 
ning from July 2000 to February 2006, in which XMM J 1046 
was visible. The event files were processed applying standard 
procedures with SAS 7.1 .20. The MOS and pn observations were 
reduced using emchain and epchain, respectively, applying de- 
fault corrections. The event lists were filtered for intervals of 
high background activity as well as to retain the pre-defined pat- 
terns corresponding to single, double, triple and quadruple pixel 
events, for the MOS observations, and single and double pixel 
events, for the pn observations, as these have the best energy 
calibration. Source and background events were extracted us- 
ing circular regions of radii ~ 25 arcsec (centered on the position 
of the X-ray source) and ~ 50 arcsec, respectively. Background 
regions were defined on an area free of sources in the same 
CCD and roughly at the same distance from the readout node as 
the source region. We restricted our analysis to the 0.15-3keV 
energy range. Whenever possible for a given observation, data 
from all EPIC cameras were analysed simultaneously to better 
constrain the spectral parameters. 



^ http : //xmmssc-www. star . le . ac . uk/Catalogue/xcat_ 
public.2XMMp.html 



^ http://xmin.esac.esa.int/sas 
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Table 2. X-ray spectral analysis of INS candidate XMM J1046. Only the set of 10 best XMM-Newton observations (see text) 
and the Chandra one are shown. The best fit parameters (derived using the maximum likelihood C statistic) listed are for an 
absorbed blackbody model for which all parameters were allowed to vary freely (middle column). The right column shows the 
blackbody temperatures obtained when the column density is held constant to the mean value, A^h = 3.5 x 10^' cm"^. Errors are 
1 cr. The observed flux fx,i4 refers to ranges 0.15-3keV and 0.5-3 keV for XMM-Newton and Chandra data, respectively, in units 
of 10 erg s ' cm"^. The "goodness-of-fit" is derived from a number of 1000 Monte Carlo simulated spectra. 



The Chandra observation was analysed using CIAO 4.0.1 
and CALDB 3.4. We processed the event files with the task 
acis_process_events, also applying default corrections. Since 
the observation was taken in VFAINT mode, we cleaned the 
ACIS background while processing the event file. We checked 
for the presence of known processing offsets using the aspect 
calculator tool available on the CIAO web pages. Finally, we se- 
lected events corresponding to grades 0, 2, 3, 4 and 6 (in ASCA 
terminology) and applied the good time intervals, also filtering 
for periods of background flares. For the spectral analysis, we 
used the 0.5-3 keV energy band due to the molecular contam- 
ination that is degrading the quantum efficiency of the ACIS 
front-illuminated chips, at energies below 0.5 Source and 
background regions were extracted applying similar criteria as 
for XMM-Newton data. In Table [1] we list, for every observation, 
the effective exposure times, off-axis angles, the number of ex- 
tracted source photons for the spectral analysis and the source 
count rate, corrected for vignetting, in the given energy bands. 



2.2. Spectral analysis 

EPIC and Chandra images do not reveal any particular back- 
ground enhancement close to the X-ray source. However, the 
Chandra data, which benefit from a lower instrumental back- 
ground, do show some weak diffuse extended emission (a fila- 
mentary structure roughly consistent with some large-scale neb- 
ulosity seen in the opticafl). The brightest part is located 5 to 
7 arcmin south-east of the X-ray source. This filament is not de- 
tected on the XMM-Newton images, although it is marginally in 
the field of view of the EPIC cameras. None of the background 
regions used in the spectral analysis overlaps with this diffuse 
X-ray emission. Therefore, it is expected that the noise for the 
EPIC data is dominated by the (position independent) instru- 
mental background. The noticeable vignetting and spreading of 
the PSF prevailing at the large off-axis angles where the source 
stands in many XMM-Newton observations both contribute to the 



* http : // cxc . harvard . edu/ciao4 . 8/index . html 
' http : //cxc . harvard . edu/ proposer/POG 
ESQ MAMA-R digitized plate. This filament is not in the field of 
view of our ESO/SOAR images. 



difficulty in measuring the spectrum of XMM J 1046 in these 
cases. 

Spectra were binned requiring a different minimum num- 
ber of counts per energy bin, depending on the total num- 
ber of source counts. Spectra extracted from the shortest expo- 
sures have at least 5 counts per bin. Using XSPEC 12.41 we 
tested different models (blackbody, power law, bremsstrahlung, 
Raymond-Smith, . . . ), allowing the fit parameters to vary freely. 
Due to the low nu mber of coun ts, we applied the maximum like- 
lihood C statistic (ICashll 19791) in order to derive the best fit pa- 
rameters and their uncertainties. The quality of each fit ("good- 
ness") is estimated by means of Monte Carlo simulated spec- 
tra, drawn from the best fit model, and the distribution of corre- 
sponding C fit statistics. If most (> 50%) of the simulated spec- 
tra has a smaller fit statistic than the current model, then it is 
unlikely that the observed data were drawn from the model. 

The spectrum of XMM J 1046 is always best fitted by a sin- 
gle absorbed soft blackbody; other models (in particular, power 
law, optically thin thermal plasmas and magnetized neutron star 
atmospheres) invariably result in worse fits and there is no ev- 
idence for an additional (hard) component. In general, the X- 
ray emission of radio pulsars results f rom the sum of thermal 
and non-thermal components (see e.g. iKaspi e t al. 2004, for a 
review). A power law component usually dominates the X-ray 
emission of young pulsars (age ^10^ yr) while old pulsars (age 
S few Myr) exhibit a weak thermal component, probably orig- 
inating from small heated polar caps, in addition to a domi- 
nating power law. On the other hand, the X-ray emission of 
the middle-aged pulsars (age ~ few 10^ yr) known as "Three 
Musketeers" is clearly dominated by soft blackbody components 
supe rposed on a high en ergy tail with photon indexes 1.7-2.1 
(e.g. 'De Luca_eLaL||2005'). The addition of a power law compo- 
nent with similar photon indexes to the blackbody model does 
not improve the fit of XMM J 1046. Any power law component 
contributes at most 4% (3 o" confidence level, 0.5-10 keV range) 
to the source luminosity. Therefore, a power law flux as low, rel- 
ative to the thermal component, as that of the Three Musketeers 
(0.3-1.7%) would not be detectable in XMM J1046 with the 
available data and can not be presently excluded. However, an X- 



http : //heasarc . gsf c . nasa . gov/docs/xanadu/xspec 
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ray spectrum dominated by a non-thermal component is clearly 
ruled out. 

The M7 show rather stable spectral and timing properties. In 
particular, the X-ray brightest source RX J 1856.5-3754 shows 
very constant flux and spectral properties over several years. 
Its X-ray spectrum is remarkably well re produced by a black - 
body with no significant deviations (e.g. iBurwitz et al]|2003l) . 
On the other hand, the second brightest source RX J0720.4- 
3125 is the only one among the M7 that has been shown to un- 
dergo long-term variations in its spectral parameters, at more or 
less constant flux. This behavi our is either possib l y cyclic and 
related to the star precession (iHaberl et all 120061: iHohle et al.1 



l2008 l and references therein) or impulsive, being interpreted as 
a sudden change on the neutron star surface accompanied with 
a simultaneous torque, cau sed by an accretion or glitch episode 
dvan Kerkwiik et alj|2007l) . 

Considering the 16 XMM-Newton observations, single ab- 
sorbed blackbodies have temperatures and column densitie^ 
ranging from 76+]^ to 160+|eV and from (l.4+i i) x 10^' to 

(8.5;^2 s) ^ 10^' cm"^; errors are Icr. In some cases, the small 
number of counts prevents a well constrained fit. In order to 
check the significance of the variations of the spectral param- 
eters of XMM J 1046, we performed & x^" test on the 16 XMM- 
Newton observations, assuming that they do not change in time 
and have values equal to their respective weighted means kT - 
113 + lleVandA^H = (3.6 + 0.9)x 10^' cm^^ (hereafter, reported 
errors on the weighted means are 3 cr, unless otherwise noted). 
We found that kT and A^h are consistent with being the same in 
all observations at an acceptable but low confidence level for the 
temperature (x\j ~ 26 a\\dx\^ ~ 19 for 15 degrees of freedom, 
respectively; the probabilities of getting a larger are ~ 4% and 
~ 19%). However, we note that among the data sets with the 
largest values of x'^ there are some observations for which the 
source elongated point spread function (PSF) overlaps the CCD 
gaps of the EPIC instruments. The correction for the missing 
part of the energy-dependent PSF might not be well calibrated, 
in particular at large off-axis angles (XMM J 1046 is located at 
~ 9 arcmin, on average, but its off-axis angle can be as large as 
~ 15 arcmin, see Table [T]). This would require the inclusion of 
an additional systematic error which would have the effect of 
lowering the x^ values. On the other hand, if we only consider 
the 10 observations for which the source is not located close to a 
gaf0, we find that the temperature and column density are steady 
over the six-year time interval at rather high confidence levels: 
kT = 117+14eV(23%)andA?H = (3.5 + l.l)x lO^' cm^^ (44%), 
see Table 121 

Alternatively, if the column density is held constant to the 
mean value A^h = 3.5 x 10^' cm"^, the computed range in kT is 
narrower, 1 1 l-132eV, and the parameters are better constrained 
(1 cr errors correspond to ~ 5% of the best values against ~ 20% 
when A^H is free to vary, see Table|2]for the individual fits). In this 
case, the variations in kT are statistically significant (99.26%). 
We note, however, that the sample of analysed observations is 
highly heterogeneous and thus subjected to systematic uncer- 
tainties. In particular, the different observing conditions - dis- 
tance to the optical axis, observing modes and cameras - as well 
as calibration uncertainties at energies below 0.5 keV especially 




Fig. 1. Contour plot corresponding to the overall X-ray fit of 
source XMM J 1046, using the 10 best XMM-Newton observa- 
tions (kT = 122+1° eV ^jj^ ^ 3 g+0.6 ^ jq2i ^^^^-2. ^^^^^ 

3 cr). The kT -Nh parameter space is shown for 1 cr, 2 cr and 3 cr 
confidence levels. The individual best fits of Table |2] are shown 
as crosses with 1 cr eiTor bars. 



for the MOS camera£3> do not allow one to draw definite con- 
clusions on the source intrinsic spectral variability. Hereafter, we 
adopt the weighted mean blackbody temperature and column 
density of the 10 best XMM-Newton observations for the pur- 
pose of further discussion on the source properties. In Fig.[T]the 
individual best fits are plotted together with the A^h xkT contours 
obtained fitting all data with one single blackbody model. 

A slightly more absorbed blackbody fit is found for the 
Chandra observation (Table |2]i. We note, however, that the re- 
sponse of the front-illuminated ACIS chips is not best suited to 
derive the spectral parameters of a source as soft as XMM J1046. 
This, and the fact that only photons with energies above 0.5 keV 
were considered, act together to cast some doubts on the column 
density derived in the blackbody fit. In Fig.|2]we show the spec- 
tral fits of the two data sets with the best signal-to-noise ratio 
among the analysed data. 

The weighted mean of the 0.15-3keV observed flux is 
fx = (9.7 + 0.5) X 10"'4 ej-g g-i cm^^, considering the 16 XMM- 
Newton observations; the computed range is (7.7 + 0.7) x 10"'^ 
to (l2.9;^[ g) X lO"'"* erg s"' cm"^ (1 cr). Although a constant flux 

is not statistically acceptable (x^ = 43 for 15 degrees of free- 
dom), a long term variability is excluded once the fluxes in the 
EPIC cameras are considered separately: the averaged flux in 
the pn camera is systematically fainter than that measured in the 
MOS cameras but both are consistent with a constant flux (at 
confidence levels of ~ 19% and ~92%, respectively). The same 
~ 13-20% discrepancy between the two instruments is seen in 
the three observations where both cameras were simultaneously 



^ The errors on A'h reported in lPires & MotchI ( [2001 Table 2) should 
be multiplied by a factor 10. 
' Observations flagged with "no" in Table[T] 



http : //xmm2 . esac . esa . int/docs/documents/CAL-TN- 
8818.pdf 
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Fig. 2. X-ray spectra and best fits of source XMM J1046 as observed with XMM-Newton (left, OBSE) 311990101, pn and MOS2 
cameras; the two lower data points consist of two subsequent M0S2 exposures) and Chandra (right, ACIS-I). The contour plots in 
the insets show the kT - Nh parameter space for 1 cr, 2 cr and 3 cr confidence levels. 



on. According to recent cross-calibration studie^Hl the MOS 
cameras register 7-9 % higher flux than pn below 4.5 keVPI and 
this excess increases with off-axis angle. The larger discrepancy 
can thus be explained by the large off-axis angle of XMM J1046 
in the set of analysed observations. The probabilities that tem- 
perature and column density are constant also increase when the 
cameras are considered individually (a consequence of the larger 
errors), but there is no systematic relation as for the flux (i.e. the 
weighted mean values are roughly the same). 

Once again, accounting only for the set of 10 observations 
not close to a gap, the overall picture further argues against sig- 
nificant flux variations: the MOS and pn fluxes can be consid- 
ered to be constant either when the cameras are analysed to- 
gether or separately (Fig. [3j. The weighted means and respec- 
tive confidence levels for the two cameras are /x,pn - (9.4 + 
1.1) X 10-14 ergs-' cm-2 (44%) and /x,mos = (1-07 + 0.08) x 
10"'-' ergs-' cm-^ (78%). Despite of the discrepancy between 
the two instruments, the fluxes are compatible within errors. 
We adopt here the mean value of the flux as measured by the 
two cameras, which is also consistent with the value inferred 
when the cameras are analysed together fx = (1.03 + 0.06) x 
lO-'-' erg S-' cm-^. A direct comparison of the source count rate, 
corrected for vignetting, for the MOS cameras, is also shown in 
Fig. [3] 

2.3. Timing analysis 

One of the most intriguing and discernible features of the M7 
is that, when compared to radio pulsars, the neutron star spin 
periods are longer and distributed in a much narrower range. 
Six of the sources show sinusoidal X-ray pulsations with pulsed 
fractions between ~ 1% and 18%. The detection of pulsations in 
the X-ray emission of source XMM J 1046 would represent fur- 
ther evidence to unveil its nature. However, despite of the large 
number of archival observations that detected the source, only 
one (OBSID 311990101, pn) is really suitable to conduct tim- 

" S. Mateos talk in the XMM-Newton Survey Science Center consor- 
tium meeting 2008, http://xmm.esac.esa.int/external/xmm_ 
data_analysis/ssc_meeting/ agenda . php 

The discrepancy between the cameras is even larger (of ~ 12-13%) 
at higher energies, according to the work just cited. 



ing analysis - the others being either too short and having few 
source counts or, for the MOS observations, taken in full frame 
mode which has a too poor frame resolution (2.6 s) for timing 
purposes. Similarly, the Chandra data was carried out with a 
3.4 s time resolution imaging configuration. 

With the aim to search for pulsations, we converted the pho- 
ton arrival times of the pn observation from the local satellite to 
the solar system barycentric frame using the SAS task barycen. 
Source photons were extracted in a smaller elliptical region (rel- 
ative to the extracted region for the spectral analysis) in order to 
avoid background events. For the same reason and, since XMM 
J 1046 does not show counts above ~ 1.5 keV, we only consid- 
ered the 0.15-2keV energy rang e. We next searched for pulsa- 
tions using a (Rayleigh) test dBuccheri et al.lfT983l) . No pul- 
sations are found to a non-constraining 30% upper limit (3 cr), in 
the 0.073-lOOs period range; the dimness of the source would 
require a longer exposure in order to significantly constrain the 
upper limit on pulsations. 



2.4. Optical follow-up 

Follow-up optical observations performed in the B and R bands 
at the ESO Very Large Telescope (VLT) in February 2007 re- 
vealed no counterpart within ~ 4.3 arcsec (^ 5 cr) from the posi- 
tion of XMM J1046, a = 10:46:08.72, 6 = -59:43:06.4, as de- 
rived running the SAS task emidetect on the best XMM-Newton 
observation 311990101. The 90% confidence level error circle 
on the position, rgo = 1.33 arcsec, is given by 2.15 ^cr^ -H cr^ysp 
where cr is the nominal error as given by emidetect and cr^y^t 
is the systematic error on the detection position as provided 
by the 2XMMp catalogue, based on reliable cross-correlation 
with the USNO Bl.O optical catalogue. Similarly, we used the 
CIAO task wavdetect to determine the position and the 90% 
confidence level error circle using the Chandra data, which is 
fully consistent with the XMM-Newton results {a - 10:46:08.7, 
6 = -59:43:06.7, rgo = 1.96arcsec). 

The VLT observations consisted of short exposure (2xl50s 
in each filter) pre-imaging data which had the goal of select- 
ing possible optical candidates for spectroscopy or deeper imag- 
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Fig. 3. Top: Long term evolution of the observed flux (0.15- 
3keV) of the INS candidate XMM J1046. Only the set of 10 
XMM-Newton observations not close to a gap is shown. Errors 
are 1 cr confidence level. The flux is stable in both pn and MOS 
detectors over the six-year time interval although there is a ~ 1 3- 
20 % discrepancy between the measured fluxes in these instru- 
ments. Observations with short exposure times and for which 
the source is located at a large off'-axis angle are highlighted. 
Bottom: MOS on-axis count rates, discriminated by filter, for 
the same set of observations. 



insFl The optical images were astrometrically calibrated us- 
ing non-saturated 2MASS and GSC-2 catalogued stars and the 
GAIA softwareQ We obtained additional imaging on the source 
at the Southern Astrophysical Research (SOAR) 4.1m telescope 
in V and Ha (6xll20s and 6x1575 s) in February 2008, which 
also failed to detect any object. In order to compute upper limits 
for the optical counterpart of XMM J 1046, we defined the limit- 
ing magnitude of the data as the magnitude of the faintest simu- 
lated star still succesfuUy measured by means of PSF fitting, on 
the exact location of XMM J1046. For all filters, the synthetic 
star was no longer detected or it was rejected (while trying to fit 
the PSF model) when the signal-to-noise ratio was worse than 
2.5. At this confidence level (98.76%), the present upper limits 
are ^ 26, my ^ 25.5 and t 25. Although some structure 
in the Ha nebular emission is present close to the position of the 
X-ray source, it is not clear if it is related to XMM J 1046 since 
similar structures also appear in other parts of the nebula, see 

Fig.a 

After correcting for photoelectric absorption and interstel- 
lar extinction, the above limits imply log(Fx//^v) ^ ^.l^^y j. 
For an assumed A^h = (3 .5 + 1.1) x lO^^cm ^ Ay ~ 1.96(6) 
dPredehl & SchmittI [19951) and the unabsorbed X-ray flu?Q is 
Fx ~ (l-4!o8) 10"^2ej.gs-i(;jn-2 in the 0.1-12keV energy 
band. The high value of the X-ray-to-optical flux ratio prati- 
cally rules out any other possibility than an INS. For instance, 
late-type (M, G, K) stars and active galactic nuclei (AGN) 
show logarithmic X-ray-to-optical flux ratios ^ - 1 and within 



Details on the optical follow-up of this and other INS candidates 
are reported in a separate paper (Pires et al., in preparation). 

http : //star-www . dur . ac . uk/ ~pdraper/gaia/gaia .html 
In this work we adopted the notation Fx and Fy to denote the un- 
absorbed fluxes while fx is used for the observed flux. 




o 




Fig. 4. Optical Ha image of the field of XMM J 1046. An in- 
verted colour map is used, i.e. brighter objects are darker. North 
and east point upwards and to the left, respectively. The im- 
age was smoothed using a Gaussian function (0.15 arcsec). This 
SOAR observation is affected by fringing. 



-1 and 1 (e.g. iBarcons et al]|2007l and references therein), re- 
spectively, while cataclysmic variable (CV) systems and BL Lac 
objects, which are among t he most extreme clas ses of objects, 
have log(Fx/Fv) ^ 2 (e.g. Schwope et alj |l999l) . Obscured T- 
Tauri stars would exhibit much more absorbed X-ray spectra. 
We also note that the high total Galactic extinction in the di- 
rection of XMM J1046 (E(B - V) ~ 12: ISchlegel e"tanil998l) 
rules out a background AGN. Moreover, the Parkes-MIT-NRAO 
Multibeam Survey of the southern hemisphere shows no radio 
source located at ^6 arcmin f rom the position of XMM J1046 
to a flux limit of ~32mJy dWright etal.1 11994 . The Parkes 
Multibeam Pulsar Sur vey h as a sen sitivity of abotit 0.14 mJy 
for a canonical pulsar (|l^^e 2008; M anchester et aljr2005l) and 
also failed to detect the source - the nearest entry in the ATNF 
Pulsar Catalogue is AXP IE 1048.1-5937, at -32 arcmin from 
XMM J 1046, and the young radio pulsar PSR J 1052-5954, at 
~ 50 arcmin. 

This, together with its thermal, soft energy distribution, 
makes XMM J 1046 a very promising INS candidate, with over- 
all properties similar to the M7. Deep radio searches and dedi- 
cated X-ray and optical observations should be carried out in or- 
der to confirm its nature and to determine to which subgroup of 
INSs it may belong. For comparison, the identification of the M7 
RX J2143.0-I-0654 with an INS was first proposed on the basis 
of a ~ 500 net counts ROSAT PSPC spectrum a nd of follow-up 
optica l observations yielding log(Fx/F\) ^ 3 jZampieri et al.l 
2001); the nature of the source was then confirmed by a dedi- 
cated XMM-Newton observation tZane et al.ii2005i) . 

3. Discussion 

On the basis of the present data, the identification of XMM 
J1046 with an INS appears convincing. If proven correct, this 
would be the first example of a presumably radio-quiet and X- 
ray dim INS, located at a significantly greater distance than the 
M7. 
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Most probably, all cooling INSs contained in the ROSAT 
Bright S ource Catalogue have alre ady been identified 
jPopov e t al. 2000; Rutledg e et alj l2003l) . so that new sources 
must be looked for at lower fluxes. The search for "blank 
field" source in the entire ROSAT All-Sky Survey and in 
HRI pointings (e lAgiieros et aLll2006l; IChieregato et al.ll2005l; 
iTreves et al.l 12007 1) produced candidates with X-ray- to-optical 
flux ratios of ~ 10 to 100. Further X-ray and optical investiga- 
tions of one such blank field led to the discovery of Calvera, a 
likely INS (which exact natur e is, however, still un clear) with 
a very large Fx/Fy > 8700 (iRutledge et al.ll2008l) . However, 
this source, like most of the other INS candidates discovered in 
the ROSAT data at faint fluxes, is at high Galactic latitude while 
it is expected that both accreting and cooling INSs at greater 
distances should be more abundant close to the Ga lactic plane. 

Detailed population synthesis calculations dPosselt et al.l 

12008 ) show that, in general, cooling INSs at lower soft X-ray 
fluxes are also expected to be hotter than the M7 (an observa- 
tional bias due to the higher photoelectric absorptions) and still 
at small ang ular distan ce from their birth star forming regions. 
Recently, M uno et all (2008) used ~ 1000 XMM-Newton and 
Chandra archival observations covering significant part of the 
sky close to the Galactic plane {\b\ < 5°) to present constraints 
on the number of magnetar candidates. The search was sensitive 
to sources with luminosities Lx > 3 x lO^-'ergs"' and pulsed 
fractions pf > 15 %, at distances within a few kiloparsecs (see 
their Fig. 6). Interestingly, their results can also be used to con- 
strain the number of cooling INSs, since no new pulsating neu- 
tron star candidate was found in spin period range ~5 to 20 s. 
Rescaling luminosity down to the typical values inferred for the 
M7 (Lx ^ 3 X 1 0-' ' erg s" ' ) we obtain that the search presented by 
Muno et al. (2008) is sensitive to such objects up to a few hun- 
dred parsecs. Taking int o account the expected distribution o f 
cooHng INSs on the sky dPopov et alJl2005HPosselt et al.ll2008h . 
we obtain as a rough estimate that the limits by Muno et al. can 
be translated into 5 1000 detectable cooling neutron stars up to 
a few hundred parsecs. In comparison with other limit this does 
not provide new important constraints. Still, a similar search in 
slightly wider period range focused on the possibility of finding 
cooling INSs is welcomed. 

The spatial location of XMM J1046, close to the center of 
the Carina Nebula, and the derived value of the source col- 
umn density suggest that it may be physically associated with 
this giant H ii region. The Carina Nebula harbours a large num- 
ber of massive stars and ha s ongoing active star formation (see 
e.g. ISmith & Brooksll2007l and references therein). The source 
column density is a factor of ~ 10 higher than those typical 
of the M7 and is consistent with the one towards Eta Carinae 
(A^H ~ 3 X 10^' cm"2, Leutenegger et al. 2003). Eta Carinae's 
distance, measured with high accuracy through t he expansio n 
parallax of its circumstellar nebula, is 2.3 kpc (e.g.l Smithll2006h . 
The distance to XMM J1046 is then likely to be comparable. 

If the star surface is an isotropic blackbody emitter, then 

^ . 3.05 ( km kpc- 

d \10-'2ergs-icm-2/ llOOeV/ ^ 

which gives a radiation radius (as seen by an observer at infin- 
ity) of R^ ~ 6.1 km for XMM J1046, assuming that it is at the 
same distance as Eta Carinae. Although smaller than the canon- 
ical neutron star radius, such a value is in agreement with what 
is measured for the M7: their redshifted radiation radii, as de- 
rived from X-ray blackbody fits and distance estimates, are in 

Sources without an optical candidate in their X-ray error circles. 



the range of ~ 2 to 7 km. However, it is well known that larger 
emission radii (and softer temperatures) are inferred when the 
surface radiation is described by more realistic models which 
take into account the overall spectra of these sources - in par- 
ticular, by invoking a geometrically thin hydrogen atmosphere 
on top of the condensed neutron star surface (e.g. iMotch et al.l 
120031: IZane e"tani2004 IHo et al.ll2007l) . Moreover, the presence 
of small hot regions on the surface is understood in terms of 
the anisotropic heat transport that occurs in the crust of cool- 
ing neutron stars endowed with a strong toroidal magn etic field 
component dPerez-Azorm et aL l l2006l:lPage etani2007h . Recent 
investigations of neutron star th ermal evolution which account 
for these effects dAguilera et al.i r2008) confirm that coolers can 
easily have polar caps with high temperatures and small radii. 

Alternatively, the rich star forming environment of the 
Carina Nebula brings in the intriguing possibility that XMM 
J 1046 could be a much older neutron star accreting from the 
ISM, probably born outside the nebula and whose orbit is 
presently intersecting the H ii region. We note that, in this case, 
position and velocity are not expected to be correlated. The 
gas m ass in the Carina Nebula is ~ 10^ Mq dSmith & Brooksl 
l2007h which, for a typical Hii region size of ~100pc, im- 
plies an average density of ~ 10 cm"^. In fact, as reported by 
Mizutani et al. (2002), two distinct electron density components 
are detected in a 30 pc area centered on the Carina I and II H ii 
regions: a high-density («e ~ 100-350 cm^) component and an 
extended low-density (We ^ 100 cm"^) component detectable 
over the entire mapped region. If XMM J 1046 is moving inside 
the nebula, the coiTesponding increase on the column density is 
~ 1.5xlO^**cm"^pc"', assuming a typical density of ~50cm"^. 
This means that the source could be up to ~ lOpc inside the neb- 
ula and A^H would still be compatible with the measured value 
and with the one derived for Eta Carinae. 

At 2.3 kpc, the luminosity of XMM J1046 is Lx ~ (9!'^) x 

lO^'^ergs"'. Although the estimated emission radius of XMM 
J 1046 is comparable to those of the M7, the higher blackbody 
temperature is responsible for the factor of nearly 10 higher X- 
ray luminosity. If XMM J 1046 is in the accretion phase, mass 
entrainment should then proceed at a rate M = Lx/rjc^ ~ 
5 X 10'~gs"', where 77 ~ 0.2 is the efficiency. The Bondi- 
Hoyle accretion rate for a star moving through the ISM with 
particle density n is M ~ 10"«(v'io)"'' gs"', where vio is the 
velocity of the neutron star relative to the ISM in units of 
lOkms"^ For « ~ 10 - 100 cm^, as appears likely inside the 
nebula, Bondi-Hoyle accretion can produce the required lumi- 
nosity but the star should move very slowly through the gas, 
V ~ lOkms"'. Radio pulsars are known to have ver y high spa- 
tial velocity (typically v ~ 400 km s''; e.g. iHobbT et al. 20051 
iFaucher-Giguere & Kaspill2006l) . thus the chance to have an ob- 
ject so slow is very small. Morevover, several mec hanisms are 
known to inhibit accretion onto a magneti zed INS dBlaes et aTl 
fT99l iToropina et al.ll2003l: Ilkhsanovll2007h . 

Inter estingly, and i n agr eement with the results of simula- 
tions by iPosselt et al.l (l2008h . a higher temperature and greater 
distance respective to the M7 are al so observed in the only 
RRAT (J1819-1458, McLaug hlin et a l. 2007) detected up to now 
in X-rays. The XMM-Newton spectrum of J1819-1458 is wefl 
fitted by a blackbody with kT ~ 140 eV and a broad absorp- 
tion feature at ~ 1 keV - similar to those of the M7, which are 
usually interpreted as evidence for magnetic fields of ~ few 
10'^ - lO'^^G. Its DM distance is ~3.6kpc and the observed 
flux is /x ~ 1.5 X 10-1^ erg s-i cm^^ (0.3-5 keV). The M7 and 
this RRAT also share similar spin periods and period derivatives. 
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On the other hand, steady or transient radio emission are not de- 
tected in any of the M7 to a rather sensitive limiting flux of ~ 
few 1 ;U Jy ( Kondratiev et al . 2008 ) . In order to better understand 
the relations between these two classes of INSs, further, deeper 
pointings in X-rays and in radio are required to firmly assess the 
nature of this and other RRAT sources. 

4. Summary and conclusions 

Overall, the present observational picture of XMM J 1046 sug- 
gests that it has the right properties to be identified with a 
younger and more distant thermally emitting INS, possibly the 
eighth member of the ROSAT-discoveied Ml. The analysis of its 
X-ray emission, although based on archival data obtained with 
non-optimal configurations, reveals an intrinsically soft energy 
distribution, possibly variable on long time scales. However, un- 
certainties on the calibration accuracy of the large off-axis angles 
at which the source was observed may account for the variations 
of the derived spectral parameters. Its mean temperature, column 
density and 0.15-3keV observed flux sxe kT = 117 + 14 eV, 
A^H = (3.5 + 1.1) X 10^1 cm-2 and fx = (1.03 + 0.06) x 
10~'^ ergs ' cm"^, where we adopted the weighted means of 
the subsample of the 10 best XMM-Newton observations. The 
present optical limits imply a logarithmic X-ray-to-optical flux 
ratio greater than ~3.1, already high enough to safely exclude 
standard classes of X-ray emitters (late-type stars, AGN and 
CVs). The column density is consistent with the source being lo- 
cated within the Carina Nebula, roughly at a distance of ~ 2 kpc. 
Assuming that the X-ray emission is well described by the mean 
parameters given above, XMM J 1046 may be more luminous 
and youngish than the M7 and thus perhaps still close to its birth 
place. XMM J 1046 is unique in the sense that it is hotter than the 
seven nearby sources and may represent an evolutionary missing 
link between the different classes of magnetars, radio-transient 
and radio-quiet INSs. The accretion scenario would require an 
unlikely low-velocity neutron star, combined with high accretion 
efficiency, to account for the X-ray luminosity of XMM J1046. 
A systematic search for pulsations is crucial in order to establish 
the true nature of this X-ray source. 

Acknowledgements. The work of A.M. P. is supported by FAPESP (grant 
04/04950-4), CAPES (grant BEX7812/05-7), Brazil, and tlie Observatory of 
Strasbourg (CNRS), France. R.T. aknowledges financial support from INAF- 
ASI under contract AAE TH-058. The work of S.P is partially funded 
by INTAS through grant 6-1000014-5706; S.R is also grateful to the 
Observatory of Cagliari for hospitality when part of this work was com- 
pleted. The authours acknowledge the use of the ATNF Pulsar Catalogue 
(http://www.atnf.csiro.au/research/pulsar/psrcat). We thank the 
anonymous referee for useful comments and suggestions which helped to im- 
prove the paper 



References 

Agueros, M. A., Anderson, S. P., Margon, B., et al. 2006, AJ, 131, 1740 
Aguilera, D. N., Pons, J. A., & Miralles, J. A. 2008, ApJ, 673, L167 
Barcons, X., Carrera, E J., Ceballos, M. T, et al. 2007, A&A, 476, 1191 
Blaes, O., Warren, O., & Madau, P 1995, ApJ, 454, 370 
Buccheri, R., Bennett, K., Bignami, G. E, et al. 1983, A&A, 128, 245 
Burwitz, v., Haberl, E, Neuhauser, R., et al. 2003, A&A, 399, 1109 
Cash, W. 1979, ApJ, 228, 939 

Chieregato, M., Campana, S., Treves, A., et al. 2005, A&A, 444, 69 
De Luca, A., Caraveo, P. A., Mereghetti, S., Negroni, M., & Bignami, G. E. 2005, 
ApJ, 623, 1051 

Eaucher-Giguere, C.-A. & Kaspi, V. M. 2006, ApJ, 643, 332 
Haberl, E 2007, Ap&SS, 308, 181 

Haberl, E, Turolla, R., de Vries, C. P, et al. 2006, A&A, 451, L17 
Ho, W. C. G., Kaplan, D. L., Chang, P, van Adelsberg, M., & Potekhin, A. Y. 
2007, Ap&SS, 308, 279 



Hobbs, G., Lorimer, D. R., Lyne, A. G., & Kramer, M. 2005, MNRAS, 360, 974 
Hohle, M. M., Haberl, E, Vink, J., et al. 2008, ArXiv e-prints 
Ikhsanov, N. R. 2007, Ap&SS, 308, 137 

Kaplan, D. L. 2008, in American Institute of Physics Conference Series, Vol. 

983, 40 Yeai's of Pulsai's: Millisecond Pulsai's, Magnetars and More, ed. 

C. Bassa, Z. Wang, A. Gumming, & V. M. Kaspi, 331-339 
Kaspi, V. M., Roberts, M. S. E., & Harding, A. K. 2004, ArXiv Astrophysics 

e-prints 

Keane, E. E & Kramer, M. 2008, ArXiv e-prints 

Kondratiev, V. I., Burgay, M., Possenti, A., et al. 2008, in American Institute 
of Physics Conference Series, Vol. 983, 40 Years of Pulsars: Millisecond 
Pulsars, Magnetars and More, ed. C. Bassa, Z. Wang, A. Gumming, & V. M. 
Kaspi, 348-350 

Leutenegger, M. A., Kahn, S. M., & Ramsay, G. 2003, ApJ, 585, 1015 

Lyne, A. G. 2008, in American Institute of Physics Conference Series, Vol. 983, 

40 Years of Pulsars: Millisecond Pulsars, Magnetars and More, ed. C. Bassa, 

Z. Wang, A. Gumming, & V. M. Kaspi, 561-566 
Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993 
McLaughUn, M. A., Lyne, A. G., Lorimer, D. R., et al. 2006, Nature, 439, 817 
McLaughlin, M. A., Rea, N., Gaensler, B. M., et al. 2007, ApJ, 670, 1307 
Mizutani, M., Onaka, T, & Shibai, H. 2002, A&A, 382, 610 
Motch, C, Pires, A. M., Haberl, E, & Schwope, A. 2007, Ap&SS, 308, 217 
Motch, C, Pires, A. M., Haberl, E, Schwope, A., & Zavlin, V. E. 2008, in 

American Institute of Physics Conference Series, Vol. 983, 40 Years of 

Pulsars: Millisecond Pulsars, Magnetars and More, ed. C. Bassa, Z. Wang, 

A. Gumming, & V. M. Kaspi, 354-356 
Motch, C, ZavHn, V. E., & Haberl, E 2003, A&A, 408, 323 
Muno, M. P, Gaensler, B. M., Nechita, A., Miller, J. M., & Slane, P O. 2008, 

ApJ, 680, 639 

Ostriker, J. P, Rees, M. J., & Silk, J. 1970, Astrophys. Lett., 6, 179 
Page, D., Geppert, U., & Kuker, M. 2007, Ap&SS, 308, 403 
Perez-Azorfn, J. E, Pons, J. A., Miralles, J. A., & Miniutti, G. 2006, A&A, 459, 
175 

Pires, A. M. & Motch, C. 2008, in American Institute of Physics Conference 
Series, Vol. 983, 40 Years of Pulsars: Millisecond Pulsars, Magnetars and 
More, ed. C. Bassa, Z. Wang, A. Gumming, & V. M. Kaspi, 363-365 

Popov, S. B., Colpi, M., Prokhorov, M. E., Treves, A., & Turolla, R. 2003, A&A, 
406, 1 1 1 

Popov, S. B., Colpi, M., Treves, A., et al. 2000, ApJ, 530, 896 

Popov, S. B., Turolla, R., & Possenti, A. 2006, MNRAS, 369, L23 

Popov, S. B., Turolla, R., Prokhorov, M. E., Colpi, M., & Treves, A. 2005, 

Ap&SS, 299, 117 
Posselt, B., Popov, S. B., Haberl, E, et al. 2007, Ap&SS, 308, 171 
Posselt, B., Popov, S. B., Haberl, E, et al. 2008, A&A, 482, 617 
Predehl, P & Schmitt, J. H. M. M. 1995, A&A, 293, 889 
Rutledge, R. E., Eox, D. B., & Shevchuk, A. H. 2008, ApJ, 672, 1137 
Rutledge, R. E., Fox, D. W., Bogosavljevic, M., & Mahabal, A. 2003, ApJ, 598, 

458 

Schlegel, D. J., Finkbeiner, D. P, & Davis, M. 1998, ApJ, 500, 525 

Schwope, A. D., Hasinger, G., Schwarz, R., Haberl, F, & Schmidt, M. 1999, 

A&A, 341, L51 
Shvartsman, V. G. 1971, Soviet Astronomy, 14, 662 
Smith, N. 2006, ApJ, 644, 1151 
Smith, N. & Brooks, K. J. 2007, MNRAS, 379, 1279 

Toropina, O. D., Romanova, M. M., Toropin, Y. M., & Lovelace, R. V. E. 2003, 
ApJ, 593, 472 

Treves, A., Campana, S., Chieregato, M., et al. 2007, Ap&SS, 308, 167 
Treves, A., Turolla, R., Zane, S., & Colpi, M. 2000, PASR 112, 297 
van Kerkwijk, M. H. & Kaplan, D. L. 2007, Ap&SS, 308, 191 
van Kerkwijk, M. H., Kaplan, D. L., Pavlov, G. G., & Mori, K. 2007, ApJ, 659, 
L149 

Wright, A. E., Griffith, M. R., Burke, B. F, & Ekers, R. D. 1994, ApJS, 91, 111 

Zampieri, L., Campana, S., Turolla, R., et al. 2001, A&A, 378, L5 

Zane, S., Cropper, M., Turolla, R., et al. 2005, ApJ, 627, 397 

Zane, S., Turolla, R., & Drake, J. J. 2004, Advances in Space Research, 33, 531 



